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We investigate numerically at various fillings the ground state of the one-dimensional Hubbard 
model with correlated hopping x (Hirsch model). It is found that, for a large range of filling 
values n around half filling, and for repulsive Coulomb interaction u < Uc{x,n), phase separation 
at a nanoscale (NPS phase) between two conducting phases at different densities occurs when 
X > 2/3. The NPS phase is accompanied by the opening of a spin gap and the system behaves 
as a Luther-Emery liquid with dominant superconducting correlations. Close to half filling, an 
anomalous peak emerges in the charge structure factor related to the density of doubly occupied 
sites, which determines the size of the droplets in the NPS phase. For l/2<x<2/3a crossover to 
a homogeneous phase, still superconducting, takes place. 

PACS numbers: 05.30.Fk, 71.10.Fd, 71.10.Hf 



INTRODUCTION 



The subject of phase separation (PS) in strongly cor- 
related fermionic materials has been widely investigated 
in connection to various physical systems, ranging from 
high Tc material^ to cold fermionic atoms trapped in 
optical lattices.^ In particular, it has been noticed how 
PS often occurs close to the transition to a superconduct- 
ing (SC) behavior. There is still no full explanation of 
such an observation. Since in high-Tc materials the SC 
phase appears upon doping an antiferromagnetic insula- 
tor, in the past the attention of both experimentalists 
and theoreticians has been mainly focused on PS in the 
underdoped regime, between an undoped insulator and 
a doped metal. More recently, experimental compelling 
evidence has grown of the presence of two types of charge 
carriers^ in cuprate superconductors and the occurrence 
of nanoscopic phase separation of the two components 
has been investigated (see for instance Ref. 4 and refer- 
ences therein). 

On the side of microscopic theoretical modeling, corre- 
lated electronic materials are well described by the Hub- 
bard model and its extensions. In this context, PS mainly 
appears as segregation of an insulating phase (either the 
half-filled antiferromagnet or the immobile pairs) within 
a low-density conducting phase ^^i^i^ Very recent 1}^ it was 
shown that in some cases PS can occur as well as coex- 
istence of two conducting phases. In particular, a high- 
density conducting phase implies the presence of mobile 
pairs in the system. In this paper we explore the pos- 
sibility that phase coexistence of two phases of different 
densities in the charge degrees of freedom generates the 
appearance of SC order, when accompanied by the open- 
ing of a spin gap. In fact, a nonvanishing spin gap may in- 
duce phase separated droplets of nanoscale size in the sys- 
tems, and subsequent quasi- long-range correlations be- 
tween the mobile pairs of the different droplets, allowing 



the superconducting correlations to become dominant. 
The model we deal with is the one-dimensional Hubbard 
model with correlated hopping, the latter describing the 
interaction between charges located on bonds and on lat- 
tice sites. The model's Hamiltonian reads: 



Hp 



Hbc = - X] [1 - ^(^^^ + ^ 

(1) 

where creates a fermion with spin a = {t, j}, a de- 
noting the opposite of a, n^cr = 4cr^*cr is the a-electron 
charge, and {ij) stands for nearest-neighboring sites. The 
parameters u and x are the Coulomb repulsion and the 
bond-charge interaction amplitude, respectively, and the 
lower case symbols denote that the coefficients have been 
normalized in units of the hopping amplitude. Moreover, 
N is the number of electrons on the D-dimensional L- 
sites lattice, so that n = N/L is the average filling value 
per site. The model has been widely studied in the lit- 
erature (see for instance Ref. and references therein, 
as well as Refs. ITol -I21I). In particular, since Hbc is not 
invariant under particle-hole transform, it has been pro- 
posed in two dimensions by Hirsch motivated by a theory 
of hole superconductivity.^^ 

In D = 1, it has by now become clear that the basic fea- 
tures of the model are well captured in the weak-coupling 
limit (x <C 1) by the standard bosonization approach, 
and resemble those of the ordinary Hubbard model with 
rescaled interaction; while for x > | a quite different 
scenario emerges at half filling. ^^^^^ In particular, for not 
too large on-site Coulomb repulsion u < Uc{x) a new 
phase appears, characterized by slowly decaying SC cor- 
relations, and incommensurate fluctuations in the charge 
correlations: the phase was denoted as ICSS (incommen- 
surate singlet superconductor). The physical origin of 
such a change in the system has not yet been fully ex- 
plained, though a subsequent paper^ suggested it could 
be related to the presence of PS between conducting 
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phases with different Fermi momenta. 

Here we shah show -through a detailed density-matrix 
renormahzation group (DMRG) numerical study- that 
in fact the SC phase found at n = 1 for x > 1/2 survives 
for a wide range of filling values. Moreover, for x > 
2/3, it turns out to be characterized by nanoscale phase 
separation (NFS). 

In particular, in Sec. [Ill we first discuss some aspects 
of the model \n D = 1 and its known properties. We 
then proceed to the study of the filling dependence of the 
ground-state energy and chemical potential (Sec. IIIip . 
and of pair-pair and charge-charge correlations and spin 
gap (Sec. lIVp . Finally, we investigate in Sec. [Vl the 
charge structure factor, thus obtaining an estimate of the 
Luttinger exponent Kp that controls the decay of charge 
and pair correlators. The discussion of the results and 
some perspectives are given in Sec. |Vll 



II. HIRSCH MODEL IN L> = 1 

The bond-charge Hamiltonian ([Tj) has both spin [5'u(2)] 
and charge [u{l)] symmetries. It is not invariant under 
particle-hole transform, though the latter can be imple- 
mented to show that the range < x < 1 is in fact rep- 
resentative of the behavior of the model at any x value. 



Indeed, under the transformation Cjcr [sgn(2x— l)]-^c 
Hbc transforms (up to constants) as 



Hbc{x,u) ^\2x-l\HBc{^g^(2x-l)x,u) , (2) 

where a = ^^^^ and a = x^u. Hence, the interval < 
X <l/2 can be mapped into > x > — oo, the Hubbard 
model {x = 0) being representative of this regime, while 
the complementary range 1/2 < x < 1 can be mapped 
into 1 < X < +00 and one expects that the integrable 
case X = 1 (Ref. Il7|) is representative of such a different 
regime. In this respect, the point x = 1/2 assumes a 
special role in that it can be simultaneously mapped into 
the two limiting cases +oo and — oo. 

In D = 1, as far as X <C 1, the Hirsch model can be 
approached by means of the bosonization technique. 
It was found that Eq. ([T]) corresponds to an effective 
Hubbard model with rescaled Coulomb interaction Ueff 
and Luttinger parameter Kp given by 



u-\-8xcoskF 

Ueff , Kp 

1 — nx 



Vf 



(3) 



where kp 



and Vi? is the Fermi velocity. The re- 



sult implies that the metal-insulator transition at half 
filling still occurs at = 0. At the same time, it suggests 
that for n > 1, Ueff may actually become negative for 
sufficiently small > 0, so that a spin gap opens, and 
the system is expected to enter a Luther-Emery liquid 
(LEL) phase. Since in this case also Kp > 1, the LEL 
phase should have dominant SC correlations. Numerical 
investigations -both at T = (Refs. JA^ and J^) and at 




Figure 1: T = phase diagram of the one- dimensional Hirsch 
model at half filling (Ref. 12>). In each phase the gapped sec- 
tors (charge, Ac and/or spin As) are reported; SDW stands 
for spin-density waves, while BOW stands for bond ordered 
wave. 



T 7^ (Ref. confirmed the validity of this scenario 

for n > 1, in some cases also at larger x- values. 
On the other hand, at x = 1 the model acquires ex- 
tra symmetries and both the thermodynamics^^ and the 
T = phase diagram^^^ can be obtained. The latter turns 
out to differ from that of the Hubbard model in many 
aspects. First of all, at half filling the metal-insulator 
transition moves up to = 4. Furthermore, below such 
value of u di new phase characterized by the presence of 
pairs and off-diagonal long-range order appears. This 
phase survives for a large range of filling values around 
n = 1. Since the model in this case has no spin gap A^, 
the phase falls into the Luttinger-liquid (LL) class. 
At half filling and u > 0, successive numerical investiga- 
tions have shown that the metal-insulator transition in 
fact moves to positive Uc{x) values as soon as x > 1/2,— 
reaching smoothly the x = 1 value [uc{l) =4]. Moreover, 
As turns out to be open for any x 1 and u < Us{x). 
So that, at variance with the bosonization predictions 
at n = 1, for < Uc{x) < Us{x) the model behaves 
as a LEL, and a new phase appears, characterized by 
slowly decaying singlet-superconducting correlations and 
incommensurate modulations in the real-space charge 
correlations The phase is denoted as ICSS. Moreover, 
for Uc{x) < u < Us{x)^ di fully gapped phase with bond 
ordered wave order is observed. The qualitative phase 
diagram of the model at half filling as derived in Refs. 
11 and [12I is given in Fig. [TJ 



A. Connection with the Simon-Aligia model 

One may wonder why the value x = 1/2 plays a rele- 
vant role in the above results. In order to elucidate the 
issue, we observe that the properties of Hbc -and in 
general of those Hamiltonians in which the interaction is 
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local- are better understood when these are represented 
in terms of on-site projection operators. The latter are 
defined as X"^^ = \a)i{P\i^ where \a)i are the states al- 
lowed at a given site and a = {0, j, |, 2}, [|2) = | ||)]. 
In such a language the nonvanishing entries of the Hamil- 



tonian matrix representation are read directly as the non- 
vanishing coefficients of the projection operators. When 
rewritten in terms of these operators, Hbc turns out to 
be a subcase of the more general Hamiltonian H intro- 
duced by Simon and Aligia^ 



J 



H = 



E 



X 



tx.X. 



2crj^C72 



2crj^0d- 



(4) 



r 



in which tx = {l — 2x) and Sx = (l—x). Besides li, the be- 
havior of Hbc is determined by the strength of tx and Sx^ 
since the original parameter x appears in its matrix rep- 
resentation only through these coefficients. Importantly, 
X = 1/2 implies tx = 0^ whereas below and above such 
value tx changes sign. Unlike the overall sign of the hop- 
ping term, which can be easily adjusted by a particle- hole 
transform, the change in the sign of tx cannot be fixed. 
A negative tx induces frustration in the motion of pairs, 
which is also driven by the positive Sx term. In particu- 
lar, for l^a^l > 1 I -in our case x > 2/3- the mobility of 
the pairs should become favored in the system, at least 
for not too large u. Numerical results available so far re- 
fer instead mainly to the particle-hole symmetric choice 

f^x — '^X 



B. 



case and phase separation 



The integrable instances of H discussed in the litera- 
ture are the standard Hubbard case {tx = 1 = Sx)^ the 
infinite-^/ Hubbard model {tx = = Sx)^ the x = 1 sub- 
case of Hbc (Ref. dli) {tx = —1, Sx =0) and the model 
recently solved in Ref. |8| (arbitrary tx < 0, = 0), 
which could be representative of the regime x > 2/3. In 
this case it is found that -whenever u < Uc{x^n)- the 
ground state is characterized by PS in the range of filling 
values ni < n < rih- The two coexisting phases have dif- 
ferent particle densities ni < 1 and rih > I amounting to 
two spinless-fermions (SF) systems. For this reason we 
shall denote the model as 2SF: in the low-density phase 
the SF fluid consists of ni electrons moving in a back- 
ground of empty sites, whereas in the high-density phase 
it amounts to 2 — holes moving in a background of dou- 
bly occupied sites. Both phases are in principle conduct- 
ing, with incommensurate Fermi momenta /c^^ = irUa 
(where a = l,h). The two densities ni and can be 
expressed^ as functions of the parameters tx, and /i, 
as follows: 



ni = 1 arccos 

TT 



ffl\ 1 /U 

V 2 / ' ^ TT ^^^^^^ 



2\tx\ 



where fi is the chemical potential, which must satisfy the 
transcendental equation 



1 



rih - ni 



--\tx\s'm7rnh 

7T 



u{nh 



1) H — siuTrn/ 

TT 



(5) 



(6) 

For tx = —1 the x = 1 subcase of Hbc is recovered: 
in this case Uh = 2 — ni. 

The presence of PS at x = 1 has a simple expla- 
nation. Let us assume for simplicity u = 0. In the 
thermodynamic limit, the ground-state energy per site 
Eo{l) coincides with that of Ns SF {Ng being the num- 
ber of singly occupied sites) and is independent of the 
number {N^) of doubly occupied (empty) sites >^ 
Eo{l) = — 2/7rsin(7rA/'s/L). As a function of the filling, 
the absolute minima in the two homogeneous phases (the 
one consisting of singly occupied and empty sites only, 
and the one consisting of singly and doubly occupied sites 
only) are reached at quarter filling (n = 1/2 = rig) and 
three-quarter filling (n = 3/2, Ug = 1/2). At any n 
within the range ni = 1/2 < n < 3/2 = n/j,, one can 
use the Maxwell construction: Ug assumes the value 1/2, 
and (2n — l)/4 doubly occupied sites are added to the 
ground state at zero-energy cost, i.e., we are in presence 
of PS. On the contrary, at x = the ground-state energy 
is that of a system of n electrons with spin moving on a 
chain, Eo{0) = — 4/7rsin(7rn/2): in this case the absolute 
minimum is reached at half filling (n = 1) and no PS is 
observed. For continuity argument, one expects PS to 
take place also for x sufficiently close to 1, and at posi- 
tive, not too large, u. 

The findings of Refs. 11 and 12 at half filling and x > 1/2 
could be consistent with a scenario of PS. As an example, 
we note that within the 2SF model, the critical line for 
PS reads Uc{x^ 1) = 4x, which is quite close to the nu- 
merical transition line of Fig. [T]in the vicinity of x = 1. 
In order to verify our hypothesis, a thorough analysis 
of the regime n 7^ 1 is required. The available numer- 
ical simulations in this case^^^^^ address mainly the is- 
sues of spin gap opening and pair-pair correlations above 
half filling, which features one would expect from the 
bosonization approach at x <C 1. The focus of the next 
sections is instead on the ground-state energy, chemical 
potential, charge-charge and pair-pair correlations, in the 
whole range n 7^ 1. 
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III. GROUND-STATE ENERGY AND PHASE 
SEPARATION 

Thanks to the fact that the total number of electrons 
N is always a good quantum number, for every value 
of u and x we have computed numerically the ground- 
state energy density Eq in the range < n < 2 by a 
series of different runs varying N from to 2L in steps 
of 2, so that it is possible to fix in each run also the 
total magnetization to 0. We have used open bound- 
ary conditions and up to 768 DMRG states with three 
finite-system sweeps in order to improve the accuracy. 
The finite-size dependence on L has been investigated by 
means of preliminary runs at L = 10 and L = 50; we 
have checked that the essential quantitative features of 
the curves Eo{n) and dEo/dn (approximated with finite 
differences 5n = 2/L) depend weakly on L. The interest- 
ing feature that we observe for suitable values of u and 
X is that there is a range of filling values ni < n < Uh 
where the energy density is linear as a function of n. This 
means that in the grand-canonical ensemble by selecting 
a suitable value of chemical potential one could change 
the number of particles with no energy cost. In Fig. [2] we 
report the results for the "canonical" chemical potential 
OEq I dn with L = 60 sites and = 1 at different values of 
X, while in Fig. [3] we fix x = 0.8 and vary u in the range 
— 4 again with L = 60. In both cases we observe that 
there exists a fiat region of values of n corresponding to 
an infinite charge compressibility 



U= 1 



> d^Eo 



The size of this region increases with increasing x and 
diminishing while for instance at x = 0.5 and u = 
1 it is absent. For reasons which will become clear in 
the next sections, we denote such a region as SC-LEL. 
The numerical phase diagrams in the planes (n,x) at 
u = 1 and (n, u) dit x = 0.8 are plotted in the insets of 
Figs. [2] and [3] and are obtained by tracing the (discrete) 
values of n at which dEo/dn starts to be flat. We should 
observe that the energy density in the SC-LEL region 
does not show a convex shape so we may safely extract 
the transition lines from the end points of the plateaux.^ 
On the one hand, a finer inspection reveals that, due 
to numerical errors and/or finite-size effects, the edges of 
the plateaux show some rounding in certain cases. On the 
other hand, we can compare our numerical results with 
the transition lines obtained analytically in the exactly 
solvable model of Ref. 8, where = 0, by taking here 
tx = 1 — 2x, and 



dEo 
dn 



— 2cos(7rn) for n < ni 

+ 2(2x — 1) cos(7rn) for n > nh , 



(8) 



where ni and nh together with the constant value of the 
chemical potential fi in the region ni < n < n^ are found 
by solving Eqs. ([5]) and (|6|). The comparison between 
DMRG and analytical estimates of ni and nh in Fig. [3] 
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Figure 2: (Color online) Discrete derivative [Eo{n + 6n) — 
Eo{n — 6n)]/26n of the DMRG ground- state energy Eq at 
u = 1: X increases from top to bottom with reference to 
n = 1. The continuous lines represent the exact results at 
= (Ref. 8). The numerical phase diagram is given in 
the inset: red diamonds are determined from the end points 
of the plateaux, and black dots refer to the analysis of Kp 
specified in the text. In the SC-LEL phase the shaded region 
is characterized by NFS. 



reveals a very good agreement at x = 0.8, meaning that 
such a value is already representative of the physics of 
the exactly solvable point Sx = 0. In fact, at x = 0.8 one 
has tx = —0.6 and Sx = 0.2. The agreement diminishes 
for X = 0.7 (Fig. [2]) and for x = 0.6 from our finite- 
size data (L = 60) we are unable to discern a plateau of 
nonzero width. This latter fact is in accordance with the 
observation that at x = 0.7 Sx and tx are almost of the 
same order (namely, Sx = 0.3 and tx = —0.4), though 
still —tx > Sx, whereas at x = 0.6 Sx is even greater than 

The SC-LEL regions x > Xc{u) at fixed or u < 
Uc{n) at fixed x, are the candidates for the ICSS phase 
away from the half-filling situation, where the transition 
points are indeed consistent with our previous analysis 
of the ICSS region. In order to confirm this indication 
an analysis of the spin gap and of the pair-pair, charge- 
charge correlations is required. 



IV. CORRELATIONS AND SPIN GAP VS 
FILLING 

The study of the previous section produced evidence 
of the fact that the coexistence of phases at different 
densities characterizes the system's behavior in a wide 
range of filling values ni < n < nh, ioi u < Uc{x^n). 
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n 



Figure 3: (Color online) Same as Fig. [2]but with fixed x = 0.8 
and varying u (increasing from bottom to top). In the inset 
the continuous line is the analytical estimate of the phase 
diagram according to Sec. Ill B[ while the crosses ("n+a") 
mark the transition points obtained by intersecting the ana- 
lytical curves out of the plateaux with the numerical values 
of the chemical potential within the plateau. The shaded re- 
gion is characterized by NFS. The jump at n = 1 and u = 4 
corresponds to the nonvanishing charge gap of the insulating 
phase. 

In particular, at n = 1 it seems to be present at least 
for X > 2/3 within the whole ICSS phase. Still at half 
filling it was noticedi^ that, due to the presence of an 
open spin gap, both pair-pair Pr = (cJ^cJ^Ci+riQ+rt) ~ 
(c]'|Q+rT)(c]"|Q+ri) and charge-charge Cr = {nifii-^r) — 
{ni){ni-^r) (where rii = rii^ + correlations become 
dominant for u < Uc{x^l). 

In this section we then explore the filling dependence 
of Pr, Cr, and A^. In Fig. [Hwe report Pl/s and C^/s 
(between site L/3 and site 2L/3 in an open chain) for a 
given L = 120 at u = 1, x = 0.8. Two aspects emerge 
quite clearly in the data shown. First, both charge-charge 
and pair-pair correlations keep a significant track of the 
undergoing transition to the PS regime, becoming appre- 
ciably different from zero only in the range ni < n < rih- 
Second, both ni and rih are indistinguishable, within our 
numerical precision, from the values which limit the coex- 
istence of phases described in the previous section. For 
this reason we have used the same symbols. It is also 
worth noticing that Pr and Cr have the same kind of 
dependence on n, in particular exhibiting a minimum at 
n = 1 and reaching their absolute maximum for n ~ 1.3. 

According to the theory of Luttinger liquids Pr and 
Cr are expected to become dominant in presence of 




Figure 4: (Color online) Charge-charge (black circles) and 
pair-pair (red squares) correlations between site L/3 and site 
2L/3 (open chains) versus filling n. Here x = 0.8 and u = 1. 
In the two insets we have reported the finite-size spin gaps 
computed as As = Eo(Sz = 1) — Eo(Sz = 0) (Sz being the 
total z component of the spin) for (a) n = 5/6 > m (As / 0) 
and (b) n = 9/6 > nn (As ^ 0). 



an open spin gap, since their power-law decay with r 
is determined by a different exponent in the LL case 
(A5 = = Ac) and in the LEL case (A5 7^ 0). Hence, 
we may expect that A5 7^ in the same range of filling 
values in which Pl/s and C^/s are. Accurate DMRG 
simulations show that this is in fact the case. The two 
insets of Fig. [H report the finite-size scaling of Ag for two 
significant cases: (a) < n < 1 and (b) n > n^. In both 
cases the results are in agreement with the findings for 
Pr and Cr, though they differ from standard bosoniza- 
tion predictions: as explained in Sec. [Ill within the latter 
approximation at > the spin gap is found to be closed 
for any n < 1 and open for any n > 1. 
The extension of the above analysis to values of x within 
the range 1/2 < x < 2/3 shows how the results are quali- 
tatively similar. They can be resumed in the observation 
that for u < Uc{x,n) and for ni < n < Uh the system 
becomes a LEL, characterized by an open spin gap and 
dominant pair-pair and charge-charge correlations. The 
observation can be exploited to infer the numerical values 
of ni and rih also in a region where the chemical poten- 
tial does not show evidence of phase coexistence. As an 
example, from the behavior of P^ 73 (see inset of Fig. [TQ|) 
we obtained at x = 0.6 = 1 as reported in fig. [2l 
Further information on the nature of the pairs in the LEL 
phase can be gained from the analysis of the pair struc- 
ture factor P{q) = J]]^ Pr exp (zgr) (not shown), which 
turns out to be peaked for vanishing q. This fact ensures 
that the electrons forming the pairs have opposite mo- 
menta (/c, —k). In order to establish whether the LEL 
phase displays dominant superconducting {Kp > 1) or 
charge {Kp < 1) correlations, the derivation and analysis 
of the Luttinger parameter Kp are now due. 
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Figure 5: (Color online) Static charge structure factor at half- 
filling for X = 0.8 and various u. The slope of the dashed line is 
1.6. For each curve, dotted lines are drawn in correspondence 
with q — 27Tnd. 



Figure 6: (Color online) Static charge structure factor at 
X = 0.8 and u = 1.0 for several fillings n [L — 120). The con- 
tinuous lines represent the analytical results obtained within 
the 2SF model in the low- and high-density regimes at the 
same n and tx (=0.6) values. 



CHARGE STRUCTURE FACTOR AND 



In this section we deepen our analysis of the phase di- 
agram of the model under investigation by numerically 
evaluating the Luttinger parameter Kp. For the half- 
filled case, the latter was estimated in Ref. [H by fit- 
ting the long-distance behavior of the equal time charge- 
charge and pair-pair correlation functions according to 
their asymptotic behavior, as predicted by conformal 
field theory equations^^ for a spin-gapped phase 
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There, it was shown that in the ICSS phase the dominant 
correlations are the superconducting ones, i.e., Kp > 1 
(for X 0.8 and u 1.0, Kp ^ 1.3). 

An estimate of Kp can also be extracted from the study 
of the static structure factor N{q) = e^^^Cr^ since the 
Fourier transform of the non-oscillating term of Eq. (|9|) 
gives the expression 

Kp = -N{q^{)). (11) 

Here we exploit this method to characterize the (su- 
per) conducting behavior of our system in the LEL 
regime. In our DMRG calculations we simulated an open 
chain of length L = 120, taking j = L/2 to reduce the 
finite-size effects due to the presence of the borders. 

In Fig. Owe plot N{q) at x = 0.8, n = 1 and various 
u > 0. One can see that within the ICSS phase, i.e., u < 
iXc(0.8, 1) = 3.05, the slope of the charge structure factor 
is weakly dependent on li, always implying a Kp > 1: in 
this case Kp ^ 1.6 ± 0.2. The latter is an overestimation 
of the actual value of Kp ^ 1.3.^^ This is in fact a general 



Figure 7: (Color online) Static charge structure factor N{q) 
at half- fining and u = 0.0 for several x values. Inset: Kp 
derived from N{q) as a function of x for three different values 
of u. The lines are just guides for the eyes. 



feature of the method employed since N{q) is numerically 
obtained for a finite lattice and logarithmic correction 
should be included. Nevertheless the present estimation 
is consistent with the previous result. 

Figure [5] also shows that the maximum position q^/ir 
depends on u. For a LEL we expect such a maximum to 
occur for g'* = 2k where at half filling kp = 7r/2. Here 
instead the actual value of g'* appears to be determined 
by the density of doubly occupied sites n^^, as the dot- 
ted lines reported in Fig. [5] in correspondence with the 
value q = 27rnd show. This feature suggests that in the 
ICSS phase the system behaves as an effective liquid of 
Ud hard-core bosons, at least for large enough x values. 
The presence of a marked peak at g = is reduced away 
from the half-filling case, as can be seen in Fig. [6l where 
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Figure 8: (Color online) Static charge structure factor at fill- 
ing n = 5/12 and = 1.0 plotted for several values of x. Inset: 
static charge structure factor at half filling in the insulating 
regime {u = 4.0) for several values of x. 



Figure 9: (Color online) Kp at x = 0.8 and = 1.0 as a 
function of the filling n [L — 120). The line is just a guide 
for the eyes and the crosses mark ui and Uh as determined in 
Sec. im 



one observes the different behavior of N{q) as a function 
of the filling. Remarkably, in the low- and high-density 
regimes the results coincide with those obtained analyti- 
cally for the 2SF model at the same tx value and Sx = ^ 
{Kp = 1/2 and q* /ir = ^kp)- Whereas, within the phase 
characterized by Kp > 1, the height of the maximum 
decreases and the width increases moving from half fill- 
ing, to disappear completely in proximity of ni and Uh- 
The maximum at ~ '^nnd disappears also at half fill- 
ing by sufficiently lowering x. As shown in Fig. [71 for 
0-5 ^ ^ ^ 0-'^ it correctly moves to = tt, with Kp > 1. 

The effect of the bond-charge interaction x on N{q) 
is visible also outside the LEL regime. Figure [8] shows 
how, by increasing x, there is a crossover both in the 
conducting regime at low density n = 5/12 and in the 
insulating regime at half filling with strong interaction 
u = 4 (inset). In the LL case it is seen that the dominant 
modulation in N{q) moves from the Hubbard type case 
with spin, qn = to the value in the spinless fermion 
case ggp = Akp^ Kp varying correspondingly. In the in- 
sulating case (inset), one observes the crossover between 
antiferromagnetic (x = 0, > 0) and the paramagnetic 
(x = 1, 1^ > 4) insulator: the charge structure factor 
is increasingly suppressed by enhancing the bond-charge 
parameter x. To summarize, in both cases increasing x 
drives the system toward a SF picture. 

An exhaustive study of the static charge structure fac- 
tor has been performed, so that by using Eq. (pT]) we 
are able to follow the behavior of Kp across the transi- 
tion to the SC-LEL regime both at half filling (inset of 
Fig. [7j) and with varying the filling (Fig. [9|). At n = 1 
the case = is peculiar in that the SC phase is en- 
tered from a conducting phase (with Kp = 1), while at 
any {) < u < Uc{x,n) the system is otherwise insulating 
{Kp = 0). 

The filling dependence of Kp is shown in Fig. [9] at 
X = 0.8. Remarkably, Kp has the same qualitative behav- 



ior of pair-pair and charge-charge correlations (see Fig. 
HI, in agreement with that expected from the conformal 
field theory equations ([9|) for a LEL. In particular, it is 
seen that all the three quantities reach their maximum 
for n ~ 1.3. Moreover Kp > 1 in the same region in 
which Pl/3 and C^/s are nonvanishing, the latter being 
very well approximated by the range ni < n < Uh^ with 
ni and as calculated for the 2SF model (same tx = 0.6, 

Sx=0). 

The situation is slightly different at x = 0.6. As shown 
in Fig. [TOl both Kp and Pl /s again have the same qual- 
itative behavior, so that also in this case we can safely 
estimate ni and as the values of n at which, for in- 
stance, Kp becomes smaller than 1. In this way we ob- 
tained the points reported in the phase diagram (x, n) 
also at X = 0.6 (see Fig. [2]), in which case the analysis 
of the ground-state energy was not conclusive about the 
presence of PS. Notice that the value n/j, = 2 at x = 0.6 is 
substantially greater with respect to the one obtained at 
X = 0.7 by the analysis of the chemical potential. This 
again could be a signal of the crossover to a regime in 
which the physics of the weak coupling limit (x <C 1) 
begins to emerge, merging with that of the strong cou- 
pling case X 1. In fact, < 1 is characteristic of the 
latter regime. In this case both Kp and /s reach their 
maximum value for n ^ 1.1. 

Finally, at x = 0.5 and u = 1^ previous numerical 
analysesii^i^ii^ suggest that the system is already in 
the weak coupling limit. The point in the inset of Fig. [2] 
for this value is simply obtained by imposing Kp = 1 in 
Eq. (O. 



VL DISCUSSION AND CONCLUSIONS 

The analyses of the previous paragraphs concur to 
complete the scenario of the physics described by the 
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Figure 10: (Color online) Kp vs n at x = 0.6 and u = 1. 
Inset: pair-pair correlations between site L/3 and site 2L/3 
(open chains) in same case. 



Hirsch model at > 0, arbitrary filling n, and < x < 1. 
Depending on the value of the bond-charge interaction 
X, one can distinguish three regions. The results are in 
fact representative of those relative to the more general 
Simon- Aligia Hamiltonian (|4]) and can be resumed as fol- 
lows: 

• For < X < 1/2 (ta; > 0) the system behaves as a 
Luttinger Liquid (LL) for n < 1, as a Mott insula- 
tor for n = 1, and possibly there is a transition to a 
superconducting (SC) Luther-Emery liquid (LEL) 
phase for n > 1, in agreement with bosonization 
predictions. 

• For 2/?> < X < 1 (tx ^ —Sx ^ 0) one can distin- 
guish a high- {n > Uh) and a low-density in <ni < 
1) region of spinless fermion (SF) LL. These coex- 
ist up to Uc{x^n) in the range ni < n < Uh < 2 
in a phase of nanoscale phase separation (NFS, see 
below) with SC properties. Within the NFS phase 
the LEL manifests an incommensurate modulation 
related to the number Ud of doubly occupied sites. 

• For 1/2 < X < 2/3 {-s^ < < 0) we are 
in an intermediate regime. A homogeneous LEL 
phase with dominant SC correlations appears for 
u < Uc{x^ n) in a wider range ni < n <2 {ni <1). 

Equation (|3]) for Kp at x <C 1 suggests the mechanism 
driving the transition to the SC-LEL state for n < 1: Kp 
in this case is a decreasing function of x, so that the elec- 
trons progressively tend to behave as spinless fermions 
(SF), which feature holds exactly when Kp = 1/2. At 
sufficiently low densities the change induced by increas- 
ing X is seen as a crossover in the static charge structure 
factor N{q): the cusp moves from q = 2kF to q = Akp- 
Also, at n = 1 and for large enough Coulomb repulsion 
u > 4, the change appears as a crossover from an antifer- 
romagnetic to a paramagnetic insulator. The acquired 



SF nature implies that the energy as a function of the 
filling reaches its minimum value at n < < 1, so that 
for greater filling values a regime of PS is favored. This 
possibility was already proven to work exactly at Sx = On 
(Ref. 8) (x = 1), in which case also for n > Uh the system 
behaves as a SF fluid, and the regime of PS amounts to 
the coexistence of the high- and the low-density fluids. 

The analytical results obtained in that case are in good 
quantitative agreement with some aspects of the results 
reported here, at least for x > 2/3. Nevertheless, in order 
to acquire SC properties, one central feature was miss- 
ing in the exact case: the presence of the spin degrees 
of freedom. The inclusion of the latter allows for the 
opening of the spin gap: exactly at = the ground 
state is fully degenerate with respect to the spin orien- 
tation, so that the spin-gap amplitude is vanishing. As 
soon as Sx 7^ 0, the spins rearrange at short distance, 
and the PS state further lowers its energy by opening a 
spin gap. This phenomenon induces a shorter scale in 
the size of the coexisting phases. While in absence of 
the spin degrees of freedom the low- and high-density 
SF fluids would be spatially segregated; when spin is 
considered, the short-distance relative orientation of the 
electrons spin becomes relevant. The coexisting phases 
split into droplets of nanoscale size A entering the NFS 
phase: an incommensurate modulation appears in the 
liquid related to the number of droplets. This is deter- 
mined by the number of doubly occupied (empty) sites 
Nd {Ne): explicitly, L/N^ < A < L/Nd. The effect is 
particularly evident at half filling, where Nd = and 
the charge structure factor exhibits a neat maximum at 
q = 27TNd/L for x > 2/3. The maximum moves to q = tt 
for 1/2 < X < 2/3. This fact can be interpreted as a 
crossover of the system to a homogeneous phase, where 
the spin gap is open, SC correlations are dominant, and 
there is no longer a length scale determined by the size 
of the droplets. As a further confirmation, in the same 
intermediate region the maximum in N{q) is observed at 
q = 7r(2 — n) for n > 1, consistently with the SC-LEL 
phase already present in this case at x < 1/2. The ex- 
istence of droplets for x > 2/3 is also supported by the 
observation that the correct behavior of N{q)^ which is in 
principle achieved by summing Cj^ over all the range of r 
values, can already be obtained, in its essential features, 
limiting the sum to the first A neighbors (not shown). 

One recognizes a posteriori that increasing the bond- 
charge interaction x from to 1 in Hbc amounts to 
passing, in the Simon- Aligia model of Eq. (|4]), from the 
universality class of the Hubbard model {tx_= Sx = I) 
to that of the 2SF model discussed in Ref. |8| (arbitrary 
tx ^ 0, = 0). It would be interesting to exploit this 
observation within the bosonization approach, comple- 
menting the results for the two fluids of electrons with 
opposite spin at x <C 1 (Ref. fl3l ) with a bosonization 
study of the 2SF model, in the limit of weak coupling Sx 
between the two SF fluids. 

We also notice that within the NFS phase, for u < 
Uc{x^ n) the actual value of u fixes Ud and hence the size 
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of the droplets. The analysis at = shows that the 
mechanism survives and even enriches also for < 0, 
and we expect further interesting physics to emerge also 
in relation to the context of cold fermionic atoms trapped 
in optical lattices. In this respect, it has been recently 
proven^^ that the Simon-Aligia Hamiltonian (|4]) is the 
correct candidate to describe these systems in proximity 
of a broad Feshbach resonance. 
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